Optical conductivity in A3C60 (A=K, Rb) 
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We study the optical conductivity in A3C60 (A =K, Rb). The effects of the electron-phonon 
interaction is included to lowest order in the coupling strength A. It is shown that this leads to a 
narrowing of the Drude peak by a factor 1 + A and a transfer of weight to a mid-infrared peak at 
somewhat larger energies than the phonon energy. Although this goes in the right direction, it is 
not sufficient to describe experiment. 
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I. INTRODUCTION 

The optical conductivity in A3C«or|(A=K, Rb) has an 
unusual and interesting behaviour nu The weight of the 
Drude peak is reduced by one order of magnitude relative 
to the weight for free electrons with the appropriate band 
mass. Much of the missing weight appears instead in a 
"mid-infrared" structure at about 0.06 eV. This suggests 
very strong interaction effects, e.g., electron-phonon or 
Coulomb interaction. The understanding of the optical 
absorption could therefore contribute much to the under- 
standing also of other properties of A 3 Ceo- 

A3C60 has orientational disorder, with the Cgo 
molecules taking, more or less randomly, one out of 
two preferential orientations.Q This orientational disor- 
der leads to a substantial modification of the optical 
conductivity in one-particle calculations. For an ordered 
system, the Drude peak collapses to a (5-fuaction, while 
the disorder leads to a broad Drude peakE The calcu- 
lated optical conductivity furthermore shows a structure 
at somewhat larger, energies than the experimental mid- 
infrared structured although the structure is less pro- 
nounced and at higher energy than in the experimental 
spectrum. More serious is, however, that the weight and 
width of the Drude peak are much larger than the exper- 
imental results. Although it is hard to separate the the- 
oretical results in a Drude and a mid-infrared structure, 
it may be estimated that the theoretical Drude width is 
more than a factor of ten too large. 

The strong reduction of the Drude width suggests very 
strong renormalization effects, e.g., due to the electron- 
phonon or electron-electron interactions. The fullerenes 
have phonons with an energy of about 0.D6 eV, which 
show a strong coupling to the electrons]! Since these 
phonons may transfer weight from the Drude peak to 
a mid-infrared structure, we here study the effect of 
phonons. 

We limit ourselves to calculating the electron self- 
energy to lowest order. This is sufficient if Migdal's 
theoremQ is valid. It is, however, questionable if this 



is true for the fullerenes, and we should keep in mind 
that higher order effects may be important. To obtain 
the optical conductivity we calculate the current-current 
response function. We can neglect vertex corrections,!^ 
since the electron self-energy is q-independent in our ap- 
proach. The current-current response function is then 
reduced to a product of two electron Green's function. 
We find that the electron-phonon coupling leads to a nar- 
rowing of the Drude peak by about a factor of 1/(1 + A), 
where A is the is the electron-phonon coupling constant. 
Although this goes in the right direction, it is by far not 
sufficient to explain the experimental data. 

In Sec. II we present the formalism and the model. In 
Sec. Ill we show the results and in Sec. IV multiplet ef- 
fects are briefly discussed. The results and other possible 
explanations of the narrow Drude peak are discussed in 
Sec. V. 



II. FORMALISM AND MODEL 



y0 



The optical conductivity is given b; 
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where 



w aft (q,u) = -i dte iw *<0|[7i(q,t),^(q,0)]|0). (2) 



Here j is the current operator and |0) is the ground-state. 
Below, we use a formalism where the electron self-energy 
is q-independent. It can then be showrlid that the ver- 
tex corrections in the current-current response function 
vanish for q — > 0, due to the odd parity of the current 
operator. We can then write the optical conductivity as 
a product of two Green's functions, only keeping a simple 
bubble of dressed Green's functions in the diagrammatic 
expansion of a. If we express the current operator as 
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the optical conductivity is given by 
Re a a0 = -Re V V v a *,v p , 
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where G is lie electron Green's function. This can be 
rewritten aslij 
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where j4 nm (w) = Im G nm {u> — i0 + )/7r and /(w) is the 
Fermi function. 

We consider the three t\ u orbitals of Ceo which are 
connected by hopping matrix elements t 



Here we only consider the number of electrons on a given 
site, and neglect the possible polarization of the charge 
on this Ceo molecule. Due to this assumption we ob- 
tain no terms in the current operator describing on-site 
transitions. Since the transitions between t\ u orbitals on 
the same site are forbidden, Eq. ( |Tl| ) is sufficient for our 
purposes. Imposing charge and current conservation 



q-j(q) = -e[H,p(q)], 
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we obtain 



q-j(q)=-^= £ W^( R i- R i) ( 13 ) 
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in the limit q — * 0. Here is the position of molecule 
i. We then obtain the current matrix elements 
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The orientational disascteO has been built into the ma- 
trix elements iy- mm / Jl3~il3 Deshpande et al. have used, 
a similar model for calculating the phonon self-energy.Ej 
We want to describe the coupling to the intramolecular 
five-fold degenerate H g Jahn- Teller modes. Due to the 
intramolecular character, the coupling has a local form. 
To describe the electron-phonon interaction, we use the 
Hamiltonian 
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where aj p h is the a phonon frequency, b rn annihilates a 
phonon with quantum mumber m, Vy ' are dimension- 
less coupling constantsl 17 H 18 l given by symmetry and <? is 
an overall coupling strength. The electron-phonon cou- 
pling constant A is then given by 



A = |iV(0)^, 

6 OJph 
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where N(0) is the density of states per spin at the Fermi 
energy. 

We now construct a consistent current operator, essen- 
tially following Ref. [01 We write the density p(i) at a 
site i as 
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The electron self-energy is calculated to lowest order 
in the electron-phonon interaction. 
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where Gmm and Djfy are the zero order electron and 
phonon Green's functions, respectively. The electron- 
phonon coupling is described by A^ m , which is expressed 

in terms of the coupling constants and the one- 

particle solutions. The interacting electron Green's func- 
tion is then obtained from the Dyson's equation 



GH = G (0 »H + G (0) ( 



(16) 



where a matrix notation has been used. 

We next discuss qualitatively how the optical conduc- 
tivity may change due to the electron-phonon interaction. 
If the band width is much larger than a typical phonon 
frequency, Migdal's theoremEl is valid. For states with an 
energy smaller than the phonon energy, the quasiparticle 
energy is then reduced by a factored 
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where A is the electron-phonon coupling. Furthermoce- 
the quasiparticle weight is reduced by the same factor.E2l 
For A3C60 it is very questionable if Migdal's theorem 
is valid, and interesting effects happen due to the fact 
that the band width is not much larger than the phonon 
frequencies.^ Nevertheless, we can expect to obtain some 
insight into the effect of the electron-phonon interaction 
by making the above assumptions, i.e., assuming that 
the electrons can be treated as noninteracting but with 
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weights and energies which are reduced by a factor Z. 
For u) > we then have 
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We replace e„ by Ze^ and (rc|j Q (q)|m) by 
Z(n|j a (q)|m)W, where the suffix refers to the non- 
interacting system. This leads to 
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where tr(°> is the optical conductivity without the 
electron-phonon interaction. For zero frequency a is un- 
changed, as it should, since the resistivity <r(0) is not in- 
fluenced by the electron-phonon interaction at zero tem- 
perature, considered here. We can see, however, that 
the energy scale is reduced by a factor of Z, and that 
the weight of the Drude peak is reduced correspondingly. 
For larger frequencies these considerations are of course 
too simple, since we then have to consider the whole 
Green's function including phonon satellites and not just 
the quasiparticlc. 
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FIG. 1. Optical conductivity cr(ui) for the phonon fre- 
quency uuph = 0.15 eV and for different electron-phonon 
coupling constants A. The figure illustrates how the Drude 
peak becomes narrower and how weight is transferred to a 
mid-infrared peak as A is increased. The inset shows a as a 
function of ui/Z, where Z — 1 + A. This illustrates how the 
width of the Drude peak is reduced by a factor of 1 + A due 
to the electron-phonon interaction. 



weight is transferred to a structure in the energy range 
0.2-0.4 eV. In the inset in Fig. [j] the same results are 
shown as a function of lu/Z. The curves now essentially 
fall on top of each other for small u>. This illustrates the 
result ( |l9| ) that the width of the Drude peak is reduced 
by a factor Z. Fig. || shows the results for a lower phonon 
frequency uj p h — 0.05 eV. The spectrum is similar as in 
Fig. [l], but the mid-infrared structure has moved to lower 
frequencies. 
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FIG. 2. The same as in Fig. [l]but as a function of uj/Z, 
where Z = 1/(1 + A). The similarity of the curves for small 
values of u)/Z but different values of A illustrates how the 
Drude peak width is reduced by a factor Z. 

From photoemission pfipr a free CgQ molecule!! and 
from neutron scattering^ it has been estimated that the 
strongest coupling is to the second lowest H ff mode at 
about 0.054 eV. From Raman scattering the strongest 
coupling j«as found for the lowest mode at about uj p h = 
0.033 eV.EJ The value u ph = 0.05 used in Fig. | should 
therefore be more realistic then the one in Fig. [I], and 
one might even argue for a still smaller value of uj p h- 
This would then tend to give an energy of the mid- 
infrared structure of the right order of magnitude, al- 
though it is still larger than the experimentally observed 
value 0.06 eV. The electron-phonon coupling is of the 
order A ~ 0.5 — l.O.cl The width of the Drude peak is 
then reduced by a factor of 1.5-2. This reduction goes in 
the right direction, but it is much too small to explain 
experiment. 



III. RESULTS 

In Fig. ^] we show the optical conductivity for a phonon 
frequency uj p h — 0.15 eV. Without electron-phonon cou- 
pling (A = 0) the spectrum shows a broad Drude peak. 
As A is increased, the Drude peak becomes narrower and 



IV. MULTIPLET EFFECTS 

An alternative mechanism for transferring weight from 
the Drude peak to the mid-infrared peak is provided by 
multiplet effects. Within the t\ u system, these are de- 
scribed by the exchange integral K between two t\ u or- 
bitals and the difference SU = U xx — U xv between the 
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direct Coulomb integral for equal and unequal orbitals. 
Here we use SU = 2K. The Cg^ molecule has a ground- 
state with spin 3/2 and states with the spin 1/2 at 3K 
and 5K above the ground-state. The valuei-of K has 
been estimated to be 0.05 eV,E3 and 0.024 eV.E3 The un- 
screened value has been found to be K. — 0.12 eV and 
within RPA screening K = 0.030 eV.EH The experience 
from atomic multiplets is that these are only weakly re- 
duced (~ 20%) relatively to what is predicted by the 
unscreened-. Coulomb integrals, both for free atoms and 
for solids rB We also find that to describe the multiplets 
in the h u — t\ u exciton, unscreened integrals give a split- 
ting of the right order of magnitude. Due to the lack of 
extensive experience for the large Ceo molecule, we nev- 
ertheless consider the whole range of estimates for the 
multiplet integrals below. If the lower values of these es- 
timates are used, the multiplet splitting is of the same or- 
der of magnitude as the energy of the mid-infrared struc- 
ture, and it is then interesting to study to what extent 
these effects can explain this structure. 

We have added a multiplet interaction to the Hamilto- 
nian in Eq. (^) 
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The simple Coulomb interaction 
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should also be added but is not considered here, since in 
simple treatments it does not give a contribution to the 
mid-infrared structure. 

We have estimated the self-energy to second order in 
8U and K and obtained 
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This has to be compared with the self-energy due to the 
electron-phonon energy, which is of the order 
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If we put K = 0.03 eV, W = 0.5 eV, A = 1 and w ph = 0.1 
eV, we find that £ E1 -P hon [ s about a factor of 35 larger 
than S Mult . This suggests that although the multiplet 
effects may transfer weight to the mid-infrared peak, the 
effect should be very small. If, on the other hand, we 
use a large value K = 0.15 eV for the multiplet integral, 
the second order self-energy due to the multiplet inte- 
grals becomes comparable to the electron-phonon contri- 
bution. In this case, however, the multiplet splitting is 



much larger than the energy of the mid-infrared peak. It 
therefore seems likely that the multiplet effects treated in 
second order theory cannot explain the energy and weight 
of the mid-infrared peak. We observe, however, the sec- 
ond order perturbation theory used here is not sufficient 
to describe the atomic limit, and that a better treatment 
conceivably could change the conclusions somewhat. 



V. DISCUSSION 

We have calculated the optical conductivity, including 
the effects of the lowest order self-energy diagram due to 
the electron-phonon interaction. This coupling reduces 
the width of the Drude peak and transfers weight to 
the mid-infrared structure at an energy somewhat larger 
than the phonon frequency. This leads to a mid-infrared 
structure with an energy of the right order of magni- 
tude, but a bit too large. The reduction of the width 
of the Drude peak goes in the right direction, but it is 
much too small. We observe that the self-energy was cal- 
culated under the assumption that Migdal's theorem is 
valid. Since Migdal's theorem is questionable for these 
systems, higher order corrections could modify these con- 
clusions. 

It is interesting that Liechtenstein et alBB found a 
rather narrow Drude peak (width ~ a few hundredths 
of an eV) in a one-particle calculation. As mentioned 
before, the Cgo molecules in A3C60 have primarily two 
different orientations. It has been found on theoreti- 
cal grounds that it is energetically favorable if neighbor- 
ing Cgo molecules have different ( "antiferromagnetic" ) 
orientations.E3ll3 The system can then be mapped onto 
a frustrated Ising model, for which the ground-state has 
a frustrated antiferromagnetic ordering.!!.] This ordering 
leads to. the narrowing of the Drude peak mentioned 
above.E3 Experimentally, a tendency to a short-ranged 
"antiferromagnetic" correlation has been foundj23 but 
under normal experimental conditions the samples are 
apparently cooled too fast to develop the long-ranged 
partial order assumed in Ref. |2^. It therefore does not 
seem likely that the partial ordering assumed in Ref. |2^ 
explains the narrow Drude peak in experimental samples 
used so far. 

It is interesting to ask what other effects may con- 
tribute to the explanation of the optical conductivity. We 
have illustrated that multiplet effects are unlikely to ex- 
plain the experimental results, at least if they are treated 
to lowest order. These systems have a strong coupling to 
a charge carrier plasmon atXLp-eV due to the oscillations 
of the three t\ u electrons. ElraEiJ In analogy with the cou- 
pling to the phonons, one may argue_that the plasmons 
have a coupling constant \ p i ~ 2.5.LJ Taking over the 
arguments from the electron-phonon coupling one might 
then expect a substantial narrowing from the coupling to 
the plasmons. This picture is, however, too simple, and a 
calculation of the electron self-energy in the so-called GW 
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approximation^ shows only a modest reduction of— the. 
band width. Actually, estimates of the specific heatO'L£l 
do not show an enhancement compared with the result 
obtained from band structure calculations, apart from 
the enhancement expected from an electron-phonon in- 
teraction with a A ~ 0.5 — 1. If these estimates are cor- 
rect, they suggest that many-body interactions do not 
reduce the dispersion in A3C60 (A-=j K, Rb). This is also 
consistent with the susceptibility,E3 which shows a very 
weak temperature dependence, implying that there is no 
narrow peak in the density of states. We should then 
not expect an explanation of the narrow Drude peak in 
terms of a mechanism which reduces the dispersion be- 
yond the reduction due to the electron-phonon interac- 
tion. Instead we should search for a mechanism which 
influences a two-particle spectrum, like the optical con- 
ductivity, without increasing the effective mass. 
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